arXiv:1501.06938vl [hep-th] 27 Jan 2015 


Dyonic Stars for Holography 

Daniel Carney 1,2 and Mohammad Edalati 1 

1 Theory Group and Texas Cosmology Center 
University of Texas at Austin, Austin TX 78712, USA 

2 Department of Physics and Astronomy and Pacific Institute of Theoretical Physics 
University of British Columbia, Vancouver BC V6T 1Z1, Canada 

We construct back-reacted asymptotically AdS 4 backgrounds with both electric and magnetic 
charge, at finite temperature and chemical potential. The system consists of a condensed bulk 
Dirac field hovering over a charged black brane. We give a detailed microscopic treatment of the 
bulk fermions, from which we derive the usual fluid approximations describing the condensate. In 
the context of holography, such a gravitational solution is dual to a (2+l)-dimensional CFT with a 
17(1) current, at finite density and strong coupling, subjected to a transverse magnetic field. 


CONTENTS 


I. Introduction 1 

II. Strategy and bulk equations 2 

III. The bulk Fermi sea 3 

A. Definition of the sea 3 

B. WKB wavefunctions 4 

C. WKB fermionic currents 5 

IV. Dyon star backgrounds 7 

V. Conclusions 9 

Acknowledgments 10 

A. Spin geometry 10 

B. Fermions in magnetic fields in flat spacetime 10 

1. Spectrum 10 

2. T = 0 statistical mechanics 11 

3. Equation of state, de Haas-van Alphen effect 13 

References 14 


I. INTRODUCTION 

Consider a high density system of fermions in 2 + 1 
dimensions with a global 17(1) symmetry, at strong cou¬ 
pling and subject to a transverse magnetic field B. As¬ 
suming such a system could be given a holographic de¬ 
scription, what is the dual gravitational solution? This 
paper proposes an answer to this question. 

To model such a system holographically, one is search¬ 
ing for a gravity solution in asymptotically AdS 4 space- 
time with a 17(1) gauge field and a weakly coupled bulk 


fermion satisfying suitable boundary conditions. Treat¬ 
ing the fermions in the probe limit corresponds to a 
boundary system in which the thermodynamics is con¬ 
trolled entirely by external sources. In order to make 
connections with realistic systems and to learn about the 
microscopies of the charges making up the state, one 
needs to include the stress-energy of the bulk fermions 
in the gravity solution. 

In this paper, we construct bulk solutions with these 
properties and with the fermions and gauge field con¬ 
tributing to the stress-energy at the same order, at ar¬ 
bitrary magnetic field strength. The solutions consist of 
a large number of bulk fermionic states localized to a 
finite region above a black brane carrying both electric 
and magnetic charge. We call these solutions dyon stars. 

Fermionic systems at finite density and magnetic field 
are critically important in condensed matter physics. 
The best-understood examples are at weak coupling, 
where the fermion spectrum consists of quantized cy¬ 
clotron orbits. This leads to a number of fascinating 
phenomenon, especially the de Haas-van Alphen oscil¬ 
lations in thermodynamic quantities in 1/B. Moreover, 
the ground state of the system is degenerate even at zero 
temperature. However, it is far from obvious what parts 
of this picture carry over to the strong-coupling regime. 
One can study the problem using the AdS/CFT corre- 
spondance, which requires solutions of the type presented 
in this work. 

Early work focused on probe fermions in electrically 
charged black brane backgrounds [1]. More recently, 
probes in magnetically charged black brane backgrounds 
[2, 3], including at parametrically large B / [i 2 [4], have 
been studied. Beyond the probe approximation, there are 
now solutions including backreaction of charged [5-8] and 
neutral [9] fermions, without bulk magnetic fields. Here 
we fill the obvious gap by constructing the solutions at 
arbitrary B / fi 2 < 1, including backreaction from charged 
fermions. We will work at small but finite temperature 
0 < T/n « BUi 2 . 

The dyon star solutions exhibit a rich and interesting 
phenomenology. They display oscillations in their density 
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profiles along the holographic direction as a result of the 
bulk dHvA effect. Like the purely electron stars [7, 10], 
the dyon stars are thermodynamically favored over pure 
black brane backgrounds at temperatures below a crit¬ 
ical temperature T c = T c (B,p), and the respective free 
energies of the solutions suggest a third-order phase tran¬ 
sition from the brane solutions to the stars. 

Yesterday [11] appeared, which overlaps strongly with 
the results presented here; the main difference is that we 
have given a detailed statistical-mechanical treatment of 
the bulk fluid. While this work was in progress the paper 
[ 12 ] appeared, which covers the same parameter space 
but uses a very different set of approximations, leading 
the authors to qualitatively different solutions than those 
presented here. 


II. STRATEGY AND BULK EQUATIONS 


We are looking for bulk solutions with AdS asymp¬ 
totics, significant fermionic backreaction, and a gauge 
field with both radial magnetic field and chemical poten¬ 
tial extending out to the boundary. At either zero or 
finite temperature, the easiest way to engineer the gauge 
field is to assume the existence of a black brane with both 
electric and magnetic charges deep in the bulk. Then we 
can start populating fermionic modes in this background 
by turning on a boundary chemical potential, search¬ 
ing for a self-consistent solution in which these fermions 
backreact on both the geometry and gauge field. 

Just as in an astrophysical fermionic star, the gravita¬ 
tional forces acting to condense the fermions must fight 
against degeneracy pressure, which can ultimately stabi¬ 
lize the fermions against collapse. This requires a high 
density of fermions, which in turn means that one can 
view them as an effective fluid obeying some particular 
equation of state. Then the gravitational dynamics in the 
region populated by the fermions can be described in the 
Tolman-Oppenheimer-Volkoff approximation [13, 14]. In 
turn, the equations of motion for the fermionic modes 
can be solved in a WKB approximation along the ra¬ 
dial direction [15]. In these approximations, the coupled 
Einstein-Maxwell-Dirac system can be reduced to a set 
of coupled ODEs along the holographic direction. 

Unlike a typical astrophysical star, here we have be¬ 
gun with a charged black brane sitting deep in the AdS 
bulk. We are anticipating that the fermions can con¬ 
dense above this horizon. If the black brane is only elec¬ 
trically charged, it is known that at T > 0 it is thermo¬ 
dynamically favored for the fermions to condense above 
the horizon, and at T = 0 the condensed fermions “fall 
into” the black brane and one obtains a Lifshitz-type ge¬ 
ometry in which the horizon has dissolved into an “elec¬ 
tron star”.[ 6 , 8 ] However, if the horizon carries magnetic 
charge, this is not possible: because the black brane acts 
as a magnetic monopole, the horizon cannot go away 


without providing some other source for the magnetic 
field. In [16] this was achieved by adding a dilaton, but 
for purposes of minimality here we will only consider the 
bulk fermion and gauge field, thus all of our solutions 
have a black brane horizon in the deep IR. 

In the bulk, we consider a (3 + l)-dimensional system 
consisting of the metric g , an Abelian gauge field A^, 

and a charged (Dirac) fermion ip', the underline differenti¬ 
ates these fields from their rescaled partners to be defined 
shortly. The Einstein-Maxwell-Dirac system is governed 
by 

s = /dV=j{T (fl _ 2A) _ _L £2 

- iip(0 - to), (1) 

where k 2 = 87tGn, R is the Ricci scalar, A = —3/L 2 
is the cosmological constant, F = dA is the U(l) field 
strength, q and to are, respectively, the charge and the 
mass of bulk spinor ip. Our conventions for the spin 
geometry, covariant derivatives, etc. are spelled out in 
appendix A. Throughout this paper, Greek indices //, v 
refer to bulk spacetime coordinates while hatted indices 
a, b refer to bulk tangent frame indices. As explained 
in the appendix, we will repeatedly use veilbeins to 
write quantities evaluated in the tangent frame; for ex¬ 
ample the local momentum is Pa = e^p^. 

We are interested in the regime in which the gauge 
and fermionic terms in this action backreact significantly 
on the gravitational solution. To do this in the large N 
limit, it is convenient to rescale various quantities. There 
are two dimensionless couplings in this system, which we 
denote 1 

7 =—, P = <n- ( 2 ) 

K 

Define dimensionless coordinates x M = xA/L. Now we 
define rescaled, dimensionless fields and derivatives 


R = L R, ip=[- 
\ 7 


L\ 


3/2 

J t 


, L 

Ip, TO = —TO 

7 


A u = —A 




L „ 1 iL 

D„, = —D„ = —d„ —— uj 


, — M 


4 7 


— fi,ab 


.A'" 1 ' - iA„. 


( 3 ) 


In terms of these rescaled fields and F = dA we have the 
action expressed in terms of our dimensionless variables: 

S=^J d 4 x\P^g^-(R +6) —-F 2 — ifi 2 ip(l/) — m)ip^. 

( 4 ) 


1 The parameter /3 2 here is morally equivalent to the phenomeno¬ 
logical parameter (3 used in [2, 6, 7]; since we are treating the 
fermions microscopically we can give a precise formula for the 
coefficient. 
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In writing the action this way, we have scaled out the 
large N ~ L/k behavior, and we will work in the regime 
where all the terms in the action contribute at the same 
order. We will often set L = 1 without loss of generality. 

We would like to work in a regime with a large density 
of occupied fermion states. This constrains our param¬ 
eters 7 , (3. First note that in holography, one requires 
L/k = 7 2 //3 1. Moreover, we will see that the WKB 
approximation for the fermions, or equivalently the lo¬ 
cally flat fluid approximation, requires 7 > l . 2 To ob¬ 
tain an order-one backreaction due to the fluid, we want 
/3 2 ~ 0(1). Putting these constraints together means 
that we are working with q 2 <C 1 . 

The equations of motion coming from (4) read 

/'"Y _ rpEM _i_ //t»F \ 

L nis 1 K 1 hv/ 

= (r) ( 5 ) 

T^D^ip = mtp. 

Here the gauge field energy-momentum tensor Tj ™ and 
the spinor energy-momentum tensor are given as op¬ 
erators by 

tEM _ p pp _ !„ p2 

W 2_ ( 6 ) 

faDv'fip + h.c ., 

while the spinor current operator is 

(7) 

The expectation values in these expressions are taken in a 
“bulk Fermi sea”, in which all the bulk fermionic modes 
with energy less than the chemical potential are popu¬ 
lated. Computing these expectation values requires us 
to specify the single-particle wavefunctions of each of the 
occupied fermion states. We do this in a WKB approxi¬ 
mation. We will construct the sea explictly and compute 
the expectation values of the current and the energy- 
momentum tensor operator in detail in section III. 

To solve the Einstein-Maxwell-Dirac system (5), we 
take the following ansatz for the metric and the gauge 
field: 

ds 2 = - f(z)dt 2 + ~^(dx 2 + dy 2 ) + g(z)dz 2 , 

z 2 (8) 

A = —h(z)dt — Bxdy. 

We have taken Landau gauge for the electromagnetic po¬ 
tential . 3 In our coordinates, the boundary of the space- 
time is at z —> 0, so the B field is oriented from the 


2 In this sense, here the TOV approximation is equivalent to the 
Thomas-Fermi [17, 18] approximation. 

3 In principle one could have allowed for B = B(z). However, this 
gives an x — z component to the stress tensor proportional to 
B'(z), which is inconsistent with our metric ansatz, so we must 
have B = constant. 


IR out to the boundary. As for the asymptotic bound¬ 
ary conditions, we take the metric to be asymptotically 
AclS 4 , i.e. f,g—> 3/z 2 as 0 —> 0. The gauge field satisfies 
A t —► —Pbdy and A y —> B^dyX, where these boundary 
quantities can be again rescaled and made dimensionless 


7 7 

Pbdy — B^fly — -jr^B 


(9) 


so that h(z) —> y, as z —> 0. Meanwhile, in the deep in¬ 
frared we assume the existence of a horizon which carries 
some electric and magnetic charge; we describe this in 
detail in section IV, where we explicitly find the gravita¬ 
tional backgrounds. 

Given the approximations described above, we will see 
that the stress tensor of the fermions takes the form of 
an anisotropic perfect fluid, 

{T* v ) = (p + P±)uy,u v + p±g + ( p z - p±)Xy.Xu, (10) 

where the local thermodynamic quantities p,p±,p z de¬ 
pend only on the holographic coordinate z and on the 
external thermodynamic parameters /i, B, and T. Here 
u^ is the unit vector along dt and the unit vector 
along d z . We give the detailed form of the thermody¬ 
namic functions below. This form of the stress tensor 
means that the Einstein-Maxwell equations reduce to a 
coupled set of ODEs along the 2 direction, which can be 
solved by numerical integration between the boundary 
and the horizon. This is the essential task of this paper. 


III. THE BULK FERMI SEA 


A. Definition of the sea 


As operators the bulk fermion energy-momentum ten¬ 
sor and current are given by the expressions in ( 6 ) and 
(7). To compute the expectation values of these opera¬ 
tors one should define a bulk state. We are after star-like 
gravitational solutions where the Dirac field has a large 
number of bound single-particle states, whose stress en¬ 
ergy make up the star. In such a state it is natural to 
compute the current expectation values by summing up 
the contributions o f all single-particle states with local 
energy u>(z) = w/\Jf{z) less than the local chemical po¬ 
tential h(z) = h{z)/\Jf(z). Here w is a rescaled energy 
ui = 7 _ 1 H. 

To proceed, we begin by writing the Dirac field in sec¬ 
ond quantization 

ip{x) = y^e~ an °*u a (x)q a + e lClat v a (yL)b\, 

a ( 11 ) 
1>\x) = E e~ inat vU*)b a + e^ut (x)at. 
a 

The single-particle wavefunctions u a (x) will be given ex¬ 
plicitly later. Using (11), one easily obtains the tangent 
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frame expressions 


a transparent fashion. Define the “magnetic algebra” 


(n = / 3 2 (^r>) = Kj E ( 12 ) 

v ^ sea 

where the sum over the sea means a sum over local fre¬ 
quencies 0 < u) a < h. Similarly, for the expectation value 
of the stress tensor, one obtains the following expressions 


i= ^ A+A '~^ 

= ~ i ]j ^( A ~ Af ) ■ 


(17) 


v J sea 

* J spa. 


u a ~h+ — u t ah a ab J u a + h.c. 


( 


Px - nh aab ) u <* + h - c - 


<n) = - £rf E^r'r^p, -Bx- 


4j~x,ab 
l 


(n) = - 2 7 ? E u « rtr ^ + h - c - 

V " spa 


4 7 


CT a6 Wa + 


(13) 


Here the momentum is = — *<9^ as usual, and we de¬ 
fined the rescaled quantities 


D P~i ~ x x 

u>=—, m = —, x = e„x=-, „ 

7 7 2 : (14) 

B = ef e?B = 2 2 B = 2 2 76. 


We have that [H, Hi] = 1 on functions of a;. In terms of 
these operators, the Hamiltonian is 


B = 




(18) 


where 

h.c. 

O = 



[(A-A^)a & + i{A + A^)ay 


ipz<J z . (19) 


Solving the Dirac equation directly is impossible given 
our general metric ansatz. The implicit 2 -dependence of 
all these terms coming from the veilbeins is what makes 
things difficult. However, things are tractable in a WKB 
approach along this axis. Let us take the ansatz (justified 
shortly) 


cj)(x) = X(x,z)Y(y,z)Z(z), (20) 


To proceed further, we need to find the single-particle 
wavefunctions u a (x). This means we need to work out 
a complete set of solutions of the Dirac equation. Given 
the general ansatz ( 8 ) for the metric and the gauge held, 
finding an analytic solution for the wavefunctions does 
not seem to be possible. However, analytic control over 
the solution can be obtained by performing a WKB anal¬ 
ysis of the Dirac equation along the radial direction, as 
a systematic expansion in the small parameter I /7 <C 1 . 
We now turn to this task. 


B. WKB wavefunctions 


where the 2 -dependence in X, Y is taken to be much 
slower than that of Z. In particular, we assume that 
the radial wavefunctions behave as 

Z ~ e i7S(z) (21) 

where S(z) is a complex-valued function that we will 
solve for in the following, and recall that we have 7 
1. The WKB approximation is that we take the 2 - 
dependence of everything in the Dirac equation (i.e. the 
metric and gauge held components) to be slow compared 
to the phase appearing here. This means that we can, 
to lowest order, commute 2 -derivatives past the veilbeins 
and terms involving the gauge held. To be precise, we 
assume that 


Consider the Dirac equation (5). We can eliminate 
the spin connection terms by rescaling the single-particle 
wavefunctions ifj = F<j> with F = [—(det g)g zz ]~ 1 ^ 4 = 
z j-~i/4' Using the time-independence of the metric, 
we can then rewrite the Dirac equation as a time- 
independent Schrodinger problem 


\d z h\, \d z e£\d\d z S\, |d 2 S( 2 )|« 7 |^S( 2 )| 2 . ( 22 ) 

In this approximation, we can work with the much 
simpler problem of the square of the Dirac equation. This 
is because given a solution y to the equation 

H\ = E 2 X (23) 


Hcj) = E<j), E = Cj — h, 


where the local Hamiltonian is given by 


h = -r* 



Ai) + im 


(15) then the function </> = (H + E)x will solve (15), to lowest 
order in the WKB approximation. We have that 


(16) 


H 2 = m 2 + 2 b 


1 

2 



(24) 


We would like to hnd a complete set of solutions, includ- where n = A^A and £ = —iT^T® = a z 0 1 - 2 x 2 is the 
ing the dispersion relation for u. One can write things in local spin operator. This is clearly block-diagonal and 
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z/zo 

FIG. 1. Local effective WKB potential V = gV, with V 
given in (27), for some low-lying Landau levels n — 0,..., 5 
in the bulk. The different colors denote different frequencies: 
black is uo = 0, blue is 0 < w < 1, and red is w ~ 1. This 
is evaluated in a backreacted star background of section IV, 
with 1. Here zo is the horizon. 

each 2x2 block is identical. The solution in a particular 
block can then be taken as 

X = N->p n (x+py/B)e iP «yZ(z) C±, (25) 

where is an eigenstate of cr z , and is the nth state 
of a harmonic oscillator with frequency b. 

All that is left to do is to work out the ^-dependence. 
Given (25), one finds the WKB problem for the holo¬ 
graphic direction is simply 

djZ(z) = 7 2 V(z)Z(z) (26) 

where the effective Schrodinger potential is given by 


We can read off the (local) energy spectrum. We label 
the eigenstates of H 2 , i.e. the functions (25), as x = 
Xpqn± in the notation of the previous few paragraphs. 
Here p indexes the momenta along the y-axis. These 
states have local energy 

E 2 n± =m 2 +q 2 +2b(n+^± ^). (29) 

Note that this is completely degenerate in the momentum 
p along the y-axis. Moreover, the state at Landau level 
n with spin down is degenerate with that at level n — 1 
and spin up, and we write 

Eqn = ? n 2 + q 2 + 2 bn = E qn _ = E 2 ^ n _^ + . (30) 

Including antiparticle states, each energy has fourfold de¬ 
generacy except the ground state n = 0 ,<r z = — 1 which 
is only doubly-degenerate. In flat spacetime, these de¬ 
generacies can be accounted for by an internal M = 2 
SUSY algebra [19, 20], and it would be interesting to 
understand how this works in AdS. 

Finally, we can act with H + E on the wavefunctions 
(25) to get a complete set of solutions to the Dirac equa¬ 
tion proper, (15). For the particle states we can begin 
with the multiplet 

V 1 _ (Xpqn-\ 2 _ (Xpq(n- 1)+) /oi \ 

A pqn l Q / 5 A pqn l Q I 

both members with E = E qn . Then, acting with H + E, 
we obtain the single-particle states 

^ qn = X' n Y p Z qn (32) 

with Y p = e llpv /\/l a box-normalized plane wave, the 
spin-oscillator wavefunctions 


V — V n ^±(z) — m 2 + 2b (n + — ± —^ — (w — h)~. (27) 

This is just a one-dimensional Schrodinger problem at 
zero energy. Note that U„_ = U( rl _i) + . To lowest order, 
the WKB solution is 


x \n = N x 

X qn = N x 


I (Eqn + m)x pqn— j 

\ytimXpq(n-l)+ - iqXpqn-) 
f {Eqn + r rn)Xpq(n— 1) + 

V -VtimXpqn- +iqXpq(n- 1)+. 


(33) 


Z(z) =e ± US(U ) 

r* , _ (28) 

S( z ) = J dz q(z ), q{z) = \J-V(z ). 

Here, the sign choice in the exponent is for incoming or 
outgoing radial modes, and q = ef q. 

The kind of gravitational solutions we will be inter¬ 
ested in are such that the effective potential (27) will ad¬ 
mit bound states in a finite region of z\ see for example 
figure 1. The bound states supported by this potential 
are the states that make up the star. The WKB momenta 
allowed in this compact region will form a discrete (if 
very dense) set, determined by a Bohr-Sonnnerfeld con¬ 
dition of the form Jqdz = N. Here N is an integer 
and 0 < zq < Z\ label the classical turning points of (27). 
These turning points are just the zeros of the potential 
since (26) is a zero-energy Schrodinger equation. 


and the radial wavefunctions Z qn = Z qn _ = Z q ^ n _x)+ 
given by (28). One can obtain quite similar formulae for 
a multiplet of antiparticle states, starting with the x’s in 
the lower block rather than the upper block. Working 
in a transverse box of coordinate area t 2 , one finds that 
normalizing the spin-oscillator wavefunctions 

SnmSij = J dxX^Xl (34) 


sets |A",y | 2 


2E qn ( E qn T 777.) 


-1 


C. WKB fermionic currents 


We need to work out the fermionic contribution to the 
stress energy. Continuing to regulate the computations 
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by considering a finite coordinate area H 2 on the x — y 
plane, the degeneracy of single-particle states is finite. 

When we sum over the local Fermi sea as in (12), we 
sum over the contributions from the states with local 
energy below the local chemical potential, 0 < E qn < h. 
In terms of the quantum numbers n, q , this means states 
with 


\q\ < q n := h 2 - m 2 , 


, Am 2 
n < uf '■= —— 
2b 


(35) 


where we defined 


as for the current, and using the magnetic algebra (17), 
one obtains 


2 n F <ln 

(Tit) = '7. ^7. Eqnd n \Zqn{z)\ 


n—0 q——q n 
p2 ™f Qn 

77: 


bn 


(T, i ) = /3 2 - 7f Y, E d n j7\Z qn (zr = (Tyy) (41) 

dZ/qn 


n—0 q——q n 
p2 -2 

{Tzz) = 7 ~rj E E d n 4^\Z qn {z)Y 

v/ n =Oq=_q n 


qn 


ra 2 = m 2 + 26?t, Am 2 = h 2 — m 2 . (36) 

Note that these quantities are local, and in particular are 
^-dependent. In regions where the fermion mass m > h, 
no states can be occupied, and the fermion stress tensor 
elements are all zero. With these conventions, the sum 
over the sea means 


E 


bP_ 

27r 


flF Qn 

EE E 


2 = 1,2 n—0 q——q n 


(37) 


Here the prefactor is from the degeneracy on the x — y 
plane: we need \p\ < M/2 so that the shifted oscillator 
is still within the transverse box, and there are M 2 /2n 
such states. The sum over n is doubly degenerate except 
for n = 0, so at each Landau level n the degeneracy is 
g n = £ 2 d n with 


dn = ^-{ 2-<M- (38) 

Z7T 

Let us begin with the current. Starting with (12) and 
computing directly with the wavefunctions (33), making 
use of the translational invariance in the x — y plane, we 
have for the local charge density 

p 2 "F 

d = 7) = P 2 —rj E E dn\Z qn (z)f. (39) 

* J n=0 q——q n 

As for the spacelike components, it is easy to check by 
direct computation that they all vanish : 4 

(J & ) = (JV) = {..r) = 0. (40) 

On to the stress tensor (13). In the WKB approxima¬ 
tion, the spin term in the stress tensor is suppressed by a 
factor of 7 _1 , so we drop it . 5 Using the same techniques 


We are interested in the situation in which a large num¬ 
ber of states are occupied, so that in particular we can ap¬ 
proximate the sum over the WKB momentum as an inte¬ 
gral. Moreover, we expect that in this limit, the fermions 
can be treated in a fluid approximation, and their equa¬ 
tion of state should be the same as in flat space, except 
with the external parameters evaluated in the local tan¬ 
gent frame. Comparing (39) to (B26) then suggests that 
we identify the local density of states, following [15], as 6 


9" p 2 

E yj \Zqn{z)\ 2 -> 



(42) 


Making this identification in the stress tensor and per¬ 
forming the dq integrals, one obtains 


(7) = 2/3 2 £ dr 


Qn 


n—0 

h F 


(T ii ) = 2p 2 Y / d n 

n 

np 

(Txx) = 4/3 2 £ dr, 

n 

nF 

(T BB )=Ap 2 J2dn 


hq n + ml In 


2 i„ / dn + h 


26nln 


q n + h 


m „ 


(43) 


hq n - ml In 


2 1 „ / dn + h 


At this stage, it is clear that we have reduced the 
fermion stress tensor to that of an axisymmetric perfect 
fluid (10). The local thermodynamic quantities appear¬ 
ing in ( 10 ) are given by the local stress tensor components 
given above, that is = e“e 5 Tlg. 

These expressions still involve a sum over Landau lev¬ 
els, which in our approximations will involve a summa¬ 
tion over rip ~ 0( 7 ) terms. These expressions are exact 


4 The fact that J z = 0 is a reflection of the WKB approximation, 
in which we neglected the electric field E z ~ h’ (z). 

5 Physically this means we are ignoring the backreaction of gravi¬ 
tational spin-orbit coupling. 


6 There is some ambiguity in whether the overall factor is absorbed 
into a radial wavefunction normalization or into this definition 
of the density of states. For our purposes it does not matter how 
we do it; even when we compute the boundary CFT correlators 
the overall wavefunction normalizations will drop out [21]. 
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to lowest order in 7 , but as a final approximation, we 
can expand the thermodynamics as a Poisson series in 
fiF ~ b~ 1 (see appendix B3 for the derivation in flat 
space). This expresses the equation of state as that of a 
free, degenerate Fermi gas at T = 0 plus contributions os¬ 
cillatory in b~ l and with amplitude supressed by a power 
of b/h 2 ~ 0( 7 -1 )- The oscillatory parts encode the fa¬ 
mous de Haas-van Alphen oscillations. Expanding, one 
obtains the zeroth-order thermodynamics 


cr 0 = 


Po = 


Po = 


/3 2 Arh 3 


37r 2 
q2 


'2 2 \ A " 4 , I h + A to \ 


8 tt 2 


24t r 2 


h(2h~ — to 2 ) Am — to 4 In 


h(2h 2 — 5m 2 ) Am + 3 In 


to j 

h + A to\ 
to / 


(44) 



FIG. 2. Local chemical potential h(z) = h(z)/yjf(z) in the 
pure black brane background (47) with T > 0. The red dashed 
line is the fixed fermion mass, black curves are the local chem¬ 
ical potential with B = 0 while blue have B > 0, and in both 
cases the horizon temperature is increasing from top to bot¬ 
tom. 


Since the equation of state at this order is independent 
of b. the local pressure po is isotropic. The leading order 
correction in the magnetic field is 


a = a 0 + 


(3 2 b 3 / 2 


00 1 

E ^3j2 Sin ( 27Ffc ^ - 


P = PO + 


^hb 3 ' 2 


fc =1 
00 


Pz=P0- 


P± = PO + 


P 2 b 5 / 2 

47 T 4 h 
p 2 b 3 / 2 
2ir 3 h 


■ E fc372 Sin ( 2nkfL F 

k—1 

v ^ 

E ^572 cos fakh F - 

k =1 


00 ^ 

A " z2 E p 72 Sin fakhp 


k =1 


(45) 


In these equations, we have a zero-temperature fluid. 
This is automatic from our definition of the bulk state 
in that we simply populated a zero-temperature bulk 
Fermi sea. However, the horizons we will assume in the 
infrared generically carry a temperature; we are work¬ 
ing in the limit in which we ignore the effect of the 
Hawking radiation on the fluid. To justify this, note 
that the local temperature felt by the fermions scales 
like Ti oca i//j,i oca i ~ 0 ( 7 -1 ), where p ioca i = 7 h is the 
unsealed, physical chemical potential, and we assume 
T//(<1 (see also [7]). 


IV. DYON STAR BACKGROUNDS 


tions (5) reduce to 

0 = - (j + - + ~)+(Pz+p)g, 

z V / 9 Z J 

0=P^ + (3 + P.-BV) 9 -i (46) 

0 = h" + ^ ( p z + p) gh! - gy/fd. 


There is by now a well-known algorithm [ 6 ] for con¬ 
structing asymptotically AdS stars satisfying these equa¬ 
tions, which we will follow here. As explained ear¬ 
lier, to get a boundary temperature and magnetic field 
we assume the existence of an electrically and magneti¬ 
cally charged planar black brane deep in the IR, say at 
z = zq > 0 . 7 In such a pure black brane background with 
T > 0, the local chemical potential h(z) is a convex func¬ 
tion equal to zero on the boundary and horizon, see fig. 
2. Thus for fixed fermion mass, one can lower the tem¬ 
perature or raise the boundary chemical potential such 
that to < h for some finite region z- < z < z+. This 
region will then have a finite density of fermion states, 
followed by another Reissner-Nordstrom solution in the 
exterior region from z = Z- out to the boundary z = 0. 
The physical parameters of the boundary theory can then 
be read off from the exterior solution. 

The interior geometry and gauge field, from the hori¬ 
zon out to the boundary of the star z_, are given by 


f(z) = 
g(z) = 


4 [1 - M 0 z 3 + (H 2 + Ql)z 4 ] 

f ,\ 4 , h(z) = po - V2QoZ, 
f(z)z 4 


(47) 

B = V2H. 


Given the expression (10) for the stress-energy of the 
fluid and the ansatz ( 8 ), one finds that the Einstein equa- 


7 These are the planar limits of the dyonic black holes discussed 
originally by Romans [22], see also [2]. 

























Here Q o and H are, respectively, the electric and mag¬ 
netic charges of the black hole, the horizon zq is de¬ 
termined by the first positive root of f(z 0 ) = 0, and 
po is a constant. Regularity of the stress tensor re¬ 
quires h(zo) = 0, that is po = V^QoZo, and the mass of 
the black hole is given by M 0 = z^ 3 [l + ( H 2 + Qq)zq] . 
Since no fermion states are occupied in this region, we 
have a = p = p z =p±=0 and one can check easily that 
(47) solves (46). 

Fix the fermion mass to, magnetic field B and bound¬ 
ary chemical potential p. Starting with sufficiently high 
T/ p, one has that h > m everywhere and the full solu¬ 
tion is simply the black brane given above. As one cools 
the system, there will be a particular value T = T c and 
radial coordinate z = z c at which 

h(z c ) = to, d z h(zc) = 0, (T = T c) (48) 

where the star is “born”. Roughly speaking, one has 
an infinitesimal shell of fermions at z = z c . Continuing 
to lower the temperature T < T c , the width of the star 
increases monotonically, occupying the region between 
the two roots z = z± of 

h{z±) = to. (T < T c ) (49) 


In the region < z < z+, a finite density of fermions 
contributes to the stress-energy and backreacts on the 
geometry. The solution in this region must be obtained 
by numerically integrating (46), subject to continuity of 
the functions /, g , h and d z h at z = z±. 

Finally, for the exterior region 0 < z < Z-, we again 
have zero fermion stress-energy, and so the metric and 
gauge fields will again be given by the Reisner-Nordstrom 
solution, although with some different parameters and an 
overall constant rescaling: 

f{z) = 4 [1 - Mz 3 + (H 2 + Q 2 )z 4 ] , 

2 c 2 (50) 

g(z ) = f( w , h{z) = c(p- V2Qz), B = V2H 
J{z)z* 


Note that we do not need to adjust the magnetic field, 
consistent with the fact that B = constant by the sym¬ 
metry of our metric ansatz. These four equations deter¬ 
mine the exterior parameters c, M, Q , p. 

The temperature as measured from the boundary is 
given by 



4-7TC 


df 


dz 


Z — 


Zo 


(51) 


The factor c is unity in the pure Reissner-Nordstrom so¬ 
lution; backreaction of the fermions sets c / 1. Given 
this equation, we can set Zq = 1 without loss of general¬ 
ity. The interior solutions are then parametrized simply 
in terms of p o and H ; by varying these parameters one 
can obtain any desired values of the ratios B/p 2 and T/p 
in the exterior. 



FIG. 3. Local charge density a(z) of some prototypical dyon 
star solutions at T > 0, with T/p increasing from top to 
bottom. One clearly sees the de Haas-van Alphen oscillations 
reflected in the stellar structure. 


In figure 3, we plot the charge density of some solutions 
of this system as a function of z. The gross features of the 
star, determined by the R-independent part of the ther¬ 
modynamics, are identical to those of the purely electron 
stars of [6, 7]. The first order correction in B/p 2 pro¬ 
vides an elegant manifestation of the de Haas-van Alphen 
effect. As the local magnetic field b = b(z) varies con¬ 
tinuously with z, the thermodynamic variables (45) are 
oscillatory in z, leading to periodic features in the stellar 
structure. 

Finally, let us study the nature of the phase transition 
between the pure black brane and the star solutions as 
the external parameters are varied. As explained in de¬ 
tail in [6, 7], the free energy of the bulk solution can be 
computed purely in terms of the boundary data, 8 

n = M- pQ- sT. (52) 

To be precise, this is the free energy per unit transverse 
area, rescaled by L/j like all the other energies consid¬ 
ered in this paper. Here s is the rescaled entropy density 
of the black brane, 

s = 27r— = 27 t —r-. (53) 

k p 

In figure 4, we show the difference in free energies be¬ 
tween the star and black brane solutions with equal exter¬ 
nal parameters. One can easily see that at low tempera¬ 
tures T < T C (B ), the star solution is thermodynamically 
favored over the pure black brane. The phase transition 
is extremely soft; in [7] it was argued to be third order, 
and the same arguments appear to apply here. 

With the chemical potential and magnetic field kept 
fixed, one can consider the zero-temperature limit. In 


The absence of a magnetic term can be understood from the 
falloffs of the gauge field: A y = Bx exactly, and there is no 
subleading behavior as a function of z. 
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FIG. 4. Top: Comparison of the free energy AQ = Q BB — 
Qstar between the star and pure black brane backgrounds. 
One sees easily that the stars are thermodynamically favored 
at sufficiently low T / fi and B// 1 2 . Bottom: Fraction of the 
bulk charge Qf due to bulk fermions. 


this limit, the black brane becomes extremal. The local 
chemical potential limits to a finite value at the hori¬ 
zon, h(zo) = fizo/y/G, and then decreases monotonically 
outward to zero at the boundary. This would suggest 
that if we take the fermion mass sufficiently small with 
fi > 0 fixed, one would have a finite density of populated 
fermions all the way down to the horizon. 

To check this idea, one can consider the fraction of the 
bulk electric charge carried by the fermions, 

Qf = Qi/J = j dz a(z). (54) 

As described in [7], in the purely electric case B = 0, as 
one lowers T//r —> 0 one finds that Qf —> 1. All the bulk 
charge is ejected from the brane into the bulk, leaving 


nothing to source the geometry behind the horizon, and 
one finds an emergent Lifshitz scaling geometry in the 
deep IR replacing the black brane. 

In our case, according to figure 4, once a finite mag¬ 
netic field is turned on, one has instead that Q / limits to 
a finite value less than one as T/fi 0. This value de¬ 
creases monotonically with increasing B/fi 2 . This means 
that the brane remains as part of the solution in the 
deep IR and continues to source both the bulk magnetic 
field and chemical potential; a finite fraction of the bulk 
charge remains behind the horizon. 


V. CONCLUSIONS 


In this paper we have proposed a gravitational dual to 
systems of fermions strongly interacting via a 17(1) cur¬ 
rent, at finite density and temperature, subjected to an 
external, transverse magnetic field. Our gravitational so¬ 
lutions, thermodynamically favored over pure black hole 
backgrounds at low temperature, consist of a condensate 
of weakly-coupled bulk fermions hovering over an elec¬ 
trically and magnetically charged black brane in the IR. 
We have called these solutions dyon stars. 

These dyon stars are solutions of Einstein gravity cou¬ 
pled only to a bulk U( 1) field weakly coupling a bulk 
fermion. This kind of minimal setup should be a robust 
approximation at low energies for a large class of field 
theories. 

We have given numerical star solutions for T > 0 and 
B > 0, and argued for a scaling limit as T —► 0 in which 
part of the bulk charge is still carried by a brane in the 
deep IR. However, explicitly constructing the T = 0 solu¬ 
tions is not possible in our setup, because the coordinate 
singularity at the horizon becomes too severe. Very near 
the horizon (for any T), various veilbeins will diverge and 
thus the WKB approximation becomes bad. One pre¬ 
sumably needs to move to some coordinates in which the 
horizon is non-singular to give the T = 0 construction. 

One particularly interesting question would be to un¬ 
derstand the relationship between the bulk “Fermi sea” 
state we have been using here, and the actual Fermi 
sea | fi) of the dual CFT constructed by occupying the 
lowest-lying CFT states. In our treatment, following 
the usual fluid approximations, we simply neglected the 
stress-energy from the parts of the fermion wavefunctions 
where the fermion is unbound. It would be interesting to 
understand such contributions in detail. 

The dynamics of strongly-coupled fermions at high 
density is of obvious and practical interest, and the addi¬ 
tion of an external magnetic field is quite relevant for real¬ 
istic lab systems. Holography continues to be a promising 
approach to such systems, and our work paves the way 
for the systematic study of CFTs with these features. 





























































































10 


ACKNOWLEDGMENTS 


We thank Willy Fischler, Sean Hartnoll, and Joe 
Polchinski for conversations, and Rob Leigh for early col¬ 
laboration. Our work at UT is supported by NSF grant 
PHY-1316033, and DC is supported by Templeton Foun¬ 
dation award 36838, by NSERC, and by the Pacific In¬ 
stitute for Theoretical Physics. 


with the Christoffel symbols; the only non-vanishing 
components in the metric (8) are 

/' 


, ,tz _ zt _ 


2 vTff’ 


(A7) 


UJ ZX = = LJ ZV = ~UJ VZ = 

x x y y 


z 2 V9 : 


and the (tangent frame) spin operator is 


Appendix A: Spin geometry 


To define the Dirac action in curved spacetime requires 
one to specify the spin geometry of the spacetime. One 
first needs to define a veilbein, which we will take to be 



V 



\ 


L^/g(z)J 


(Al) 


= 



(A8) 


For the purpose of rescaling the Dirac spinor in section 
IIIB, note that an easy calculation gives the spin con¬ 
nection’s contribution to the Dirac equation: 


rH,a^ aS = id* H-(detg)g zz }r z . (A9) 


Appendix B: Fermions in magnetic fields in flat 
spacetime 


Here fi € {t,x,y,z} is the bulk spacetime index while 
a € {t, x. y, z} is the bulk tangent frame index. As usual 
frame indices are raised and lowered with the flat metric 
7 ?afc = 9^ e & e l- We will often define “local” or “tangent 
frame” quantities by exchanging spacetime for frame in¬ 
dices; for example, the local energy and momenta are 
defined by 

Pfj. = (E,Pi) => Pa = = (£>;)• (A2) 

We define spacetime gamma matrices = e^T a as a 
section of the bulk Clifford bundle determined by some 
fixed frame gamma matrices T a . These satisfy the Clif¬ 
ford algebra 


In this appendix we review some pertinent facts about 
the dynamics of relativistic fermions in a constant back¬ 
ground magnetic field in flat spacetime. We first obtain 
the single-particle energy spectrum and eigenstates of the 
fermions. We use these to compute, in the grand canon¬ 
ical ensemble, the energy density, pressure, and charge 
density a gas of such particles. We obtain explicitly these 
thermodynamic quantities in the limit T/y, <C B/y 2 <C 1, 
up to the leading order contributions in B/y, 2 . The ze¬ 
roth order answer is that of a free degenerate Fermi gas 
while the first order corrections are periodic in 1/B. 


1. Spectrum 


{F a ,r S } = 2r] ai => {r 7 \r"} = 2 g yv . 

In this paper we will use the basis 


r* = 


i o 
0 -i 


r = 



(A3) 


(A4) 


We use hats on the Pauli matrices since they really belong 
to the tangent frame. Dirac conjugation is defined as 
usual, in terms of the spacetime gamma matrices: ip = 

We write derivatives transforming under both coordi¬ 
nate and gauge transformations, 


D n=d„- db - 


(A5) 


The second term involves the spin connection, which has 
coefficients 


Consider a Dirac fermion ip with mass m and charge 
q propagating in four-dimensional Minkowski spacetime 
with metric ?= (—1,1,1,1). The fermion is subject to 
a constant background magnetic field, which we assume 
to be in the ^-direction, F = Bdx A dy. We will work in 
Landau gauge 


A = —ydt + Bxdy 


(Bl) 


with /i the chemical potential. The fermion dynamics is 
determined from the action 


S= d 4 xiip — m)ip, 


(B2) 


where T 71 can be taken in eg. the basis (A4), and D M is 
the covariant derivative D M — iA /t . 

The energy spectrum of ip is determined by solving the 
Dirac equation 


W A1,oS — 


(A6) 


(r^ — m) ip = 0. 


(B3) 
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The problem is time-independent and we can take 
ip(x., t) = e _lart it(x). Then we can rewrite the Dirac equa- 
tion as a time-independent Schrodinger problem, 

Hu = Eu 

nr . n (B4) 

H = -r° [r(pi - Ai) + im\ , 

where as usual pi = —id-i, and we defined E = u — p. 
The spinor structure here is somewhat complicated, but 
notice that if we have a solution to the squared equation 

H 2 X = E 2 X (B5) 

then the state 

u = (H + E) X (B6) 

will solve the Dirac equation (B4). In Landau gauge 
(Bl), the square of the Hamiltonian is simply 

H 2 = m 2 + p 2 + {p y - Bx ) 2 + p\ - BY. (B7) 


potential of the cyclotron orbits exactly cancels that from 
the spin. Thus for fixed p, q 1 the energy spectrum is four¬ 
fold degenerate (counting the antiparticle states), except 
for the ground state n = 0, <r z = +1 which is only doubly- 
degenerate. 

To write the explicit solutions to (B3), it is convenient 
to write 

<BU> 

where 

O = [{A - A^)a x - i{A + A^)a y ] - iqa z . (B13) 

This operator acts on the spin-oscillator part of the wave- 
functions as 

Oip n C+ = -V2Bm/j n - iC- - iqipn C+ /T11 

, - (B14) 

OlpnC- = \J2B(n + l)^n+iC+ + iqipnC-- 


where E = — iT x F v = <r z ® 12x2 is the spin operator. This 
is clearly block-diagonal and each 2x2 block is identical. 

Solving (B5) is simple. The problem is symmetric un¬ 
der translations in both y and z, i.e. p y and p z commute 
with the Hamiltonian, so the eigenstates of H 2 are plane 
waves x ~ e lpy+zqz . Define the “magnetic algebra” 


V2B 


Px = —t] 


A + A^ + 


[A-Al], 


(B8) 


These satisfy [^4,^4^] = 1 on functions of x. The con¬ 
stant shift corresponds to a shift of the axorigin by p/B. 
Clearly the eigenstates are of the form 

X = Xp</n±( x ) = Ne l(pv+qz) tpn{x - p/B) C± (B9) 


where %j) n is the nth harmonic oscillator wavefunction, 
and is a spin- 2 : eigenstate. This is the solution in 
either of the blocks, and we can take the other block to 
simply vanish. The energies satisfy 

^ pqn± = m 2 + q 2 + 2B + - =F ^ ■ (BIO) 


The index n is known as the Landau level of the state. 

This spectrum has two notable features. First, the 
momentum along y has dropped out, so the spectrum is 
continuously degenerate in that quantum number. More 
interestingly, the spin-dependent term has exactly the 
right magnitude so that 

Eqn := + 9“ + 2 Bn = E 2 qn+ = E 2 q( ^ n _ 1 - ) _, (Bll) 

i.e. the state at Landau level n and spin up is degenerate 
with the state at level n — 1 with spin down: the dipole 


Now we can use this and the construction (B6) to get 
the full solutions. For the particle states we can begin 
with the multiplet 


Y 1 - ( Xpgn+\ 2 _ f Xpq(n-l)-\ IBIS') 

Xpqn l Q j i Xpqn l Q ) 


both members with E = E qn . Then, acting with H + E, 
we obtain the orthogonal pair of states 


/ (Bqn “t“ Xn./Xpqn+ 
yV^2 ~BnXpq{n—l)— 4” iqXpqn+ 
( E qn + m)Xpq(n- 1)- 
\/2 BnXpqn+ iqXpq(n— 1 ) 


(B16) 


One can obtain quite similar formulae for a multiplet 
of antiparticle states, starting with the x’s in the lower 
block rather than the upper block. 

From here out we will sometimes use the schematic in¬ 
dex a = {pqni} to label these states and energies. In or¬ 
der to normalize the second-quantized Hamiltonian H = 
E a ®« a a fl a, we want to normalize J d 3 xu' a up = 5 a p. 
For simplicity let us work in a box of volume V = L 3 . 9 
This sets 


|7V| 2 = [2 L 2 E qn [E qn + to)] 1 . (B17) 


2. T = 0 statistical mechanics 


Consider now a large number of these fermions. Let’s 
study the system in the grand canonical ensemble where 


9 Although we still follow the common practice of taking 
f (x)iprn (x) = Snm for the oscillator wavefunctions. 
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the chemical potential, magnetic field and temperature 
are kept fixed. 

The energy spectrum 

E 2 n = to 2 + q 2 + 2 Bn (B18) 

has a huge degeneracy, as remarked above. To regu¬ 
late the computations, note that the totally degenerate 
momentum p along the y-axis is bounded from above, 
\p\ < LB/ 2, due to the fact that the shifted center of the 
^-coordinate harmonic oscillator Xg = p/B must stay in¬ 
side the box. 10 The states with n > 1 are all doubly 
degenerate due to spin. So, at each level n, the total 
degeneracy becomes 

g n = 2irAd n , d n = (2 - (5„ j0 ) (B19) 

(2ti -y 

where we took the usual L/2it normalization for the plane 
wave factor, and the reason for defining d n this way will 
become clear shortly. 

Let us fix some notation by considering the physics 
at strictly zero temperature. For T = 0, all states with 
E qn = ojqn — g < 0 will be occupied, and none of the 
others. In each level n we define an effective mass 


Fix a chemical potential g > 0. The Fermi sea is the 
state 

im)= n 4i«>, o 824 ) 

0<cJq, </x 

where the product is ordered with energy increasing from 
right to left. Due to Fermi statistics, this state has the 
algebraic properties 

(p\ a la a '\g) = S a , a '0(c*J a < p) 

. y _D^O J 

{g\a a a\,\g) = S a ^ [1 - 0(w Q < g)\. 

When we take vacuum expecation values of bilinear op¬ 
erators like the current and stress tensor, we subtract 
the contribution from the usual vacuum state a a |0) = 0, 
and we drop any anti-particle contributions, yielding fi¬ 
nite answers. 

We are now in a position to compute the expectation 
values of the current and energy-momentum tensor op¬ 
erators. Let us start with the expectation value of the 
current. Using (B23), (B25), one easily obtains, for ex¬ 
ample, the charge density 

cr = (J 4 ) = (/#rVlM) = Y U » Ua - ( B26 ) 


to 2 = to 2 + 2Bn. (B20) 

For fixed g > 0 we have that, in a given Landau level n, 
the momentum states that are occupied are those with 

q 2 < ql ■= P 2 - m 2 n . (B21) 


Similarly, there will be a maximum Landau level with oc¬ 
cupied states, given by the floor of the continuous quan¬ 
tity 71 f- Let A m 2 = g 2 — to 2 , then let 


Uf 


Am 2 
2 B ' 


(B22) 


If the highest Landau level only has its q = 0 states occu¬ 
pied, then 71 f will be an integer equal to that level. As B 
or g vary, varies continuously and tracks the fraction 
of the uppermost Landau level that is filled. Shortly we 
will see that, in the limit where T/g < B/g 2 <C 1, the 
thermodynamic quantities are oscillatory in ~ 1/B. 

In the strict T = 0 limit, one can explicitly sum over 
the occupied states co < g to obtain the usual thermody¬ 
namic quantities. The cleanest way to do it is to compute 
the expectation of the stress energy tensor components 
in the usual Fermi sea state. Expand the Dirac field 


Here the sum over the sea means a sum over frequencies 
0 < w Q < g, i.e. 


E 


BV 


n F 

EE/ 


dq, 


(B27) 


with riF, q n defined above, and where we have taken the 
continuum limit for the 2 -axis momenta. Using the ex¬ 
plicit states u a given above, and using the fact that our 
state is translation invariant in every direction to write 
f d 3 x/V = 1, this reduces to 

n f 

cr = 2 Y dnq n - (B28) 

n—0 

Going through quite similar manipulations, one easily 
finds that (J 1 ) = 0 along the spatial axes. 

The expectation values of the stress tensor components 
are a bit more computationally involved but straightfor¬ 
ward. One obtains 



(B29) 


i>( x ) =Y e iaJ “* M a(x)a Q + e luJat v a (x)b'l, 

(B23) 

V ,t ( x ) = Y e~ lMat vi{x)b a + e Wo,t u [ a (x)a^. 


10 Another way to see this is that in axial gauge A ~ B(xdy — ydx) 
the wavefunctions are localized annuli on the plane, which must 
remain sufficiently within the box. 


where these thermodynamic functions are given by 

rip 


p 


— E dn 

n —0 
n p 

Pz = Y dn 

n —0 
n F 


p± 


= 2 Y d » 


gq n + to 2 In 
gq n - rn 2 In 
2 Bn In 


q n + M 


q n + P 

m n 


(B30) 


n —0 


q n + p 

m n 
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The sum over Landau levels is not possible to evaluate in 
closed form. However, one can work at small finite tem- 
peraute, and perform some approximations in the exact 
partition function which give tractable and interesting 
formulas. We now turn to this. 


gives the grand potential of a free Fermi gas 


fin — — 


A m 2 V 

87T 2 


2 A 1 . 2 
-pAm — - Am 
o o 


(B34) 


3. Equation of state, de Haas-van Alphen effect This leads to the zeroth-order contribution to er, p, p z 

and p±: 


In the absence of magnetic field, the energy spectrum 
is given by to 2 = in 2 + k 2 , and states whose energies 
are below the Fermi energy Wp = p 2 = m 2 + kp are 
all occupied. This implies that the B —► 0 limit should 
be taken provided that 2npB = p 2 — m 2 remains fixed. 
In this limit, the number of Landau levels with occupied 
states below the Fermi surface becomes infinite, so we can 
approximate the sum over Landau levels by an integral. 

In the limit that T/p < B/p 2 <C 1, thermodynamic 
quantities develop a part which is periodic as a function 
of 1/B. This is the de Haas-van Alphen effect. One can 
see it by thinking of the fermion states as cyclotron orbits 
on the plane, whose radii are quantized in units of the 
flux. The total thermodynamics is then the sum of a 
degenerate Fermi gas plus corrections that oscillate with 
the external magnetic field, as well as terms monotonic 
in the magnetic field. These oscillations can be easily 
seen by simply plotting the exact thermodynamics (B30) 
as a function of B. In what follows we get an analytical 
approximation for this effect. 

Working at arbtirary temperature and held, we can 
write down the exact grand potential. Taking the contin¬ 
uum limit of the 0 -axis momentum q , the grand potential 
takes the form 


Q = —TlnZ 


(B31) 


= -VTZ*n 

n—0 



1 + e ~ {u ~^ /T 


(B32) 


In order to analytically obtain the oscillatory behavior 
of the thermodynamic quantities, we can use Poisson’s 
summation formula 11 to write the grand potential f1 in 
terms of oscillating functions. Let f{n) — J dq In.... Now 
apply Poisson summation to this. 

The first term is non-oscillatory, independent of B , and 


11 This comes from writing the Dirac comb as its Fourier series, 
E^L_oo ~ n ) = oo e 27Ttkx , and integrating against a 

test function / on x £ [0, oo]. The Poisson summation identity 
is then 

- y>- 5 „,o)/(n)= / dx f(x) + 2ReJ2 / dx e 2 ” ik *f(x). 

2 n—0 J 0 fc = l J 0 




Po 

Po 


A m 3 

37T 2 


1 

87 r 2 

1 

247T 2 


m 4 In | 

' p + Am’ 

v m 

+ 3 In | 

( H + Am 

l m 


(B35) 


Since the equation of state at this order is independent 
of B , the pressure po is isotropic. 

The oscillatory integrals contain all the information 
about the magnetic held. The integrals produce both 
oscillatory and monotonic functions of B , but the latter 
are supressed by an extra factor B/p 2 , so at small B/p 2 
we get the dominant oscillatory behavior 


n 


osc 


VTB 3 / 2 

27T 2 

fc =1 


cos(27rfc?fF — 7 t/4) 
k 3 / 2 sinh(27 j 2 kpT/B) 


VB 5 / 2 
4ir 4 p 


OO 


E 


k 5/ 2 


COS 


^27rfcriF 



(B36) 


where np has been dehned in (B22), and in the second 
line we expanded around low temperatures T/p <C 1. 

Thus, at zero temperature, and to leading order in the 
magnetic held, one obtains 


B 3 / 2 AA 1 . / 7r\ 

+ -^ 3 - 2^ pp sm [2nkn F - -J 

k =1 

pB 3 1 2 1 . / 7T\ 

P=P0+ — 3 ~ 2^ Sln ( 27r ^F - J 

k =1 

B 5/2 °A i f 7r\ 

p ‘ = p °-i^2^¥K coe { 2l ' knr -4) 

A k= 1 

B 3 ! 2 2 ^A 1 . / 7 r\ 

p± = Po + 2^ Am E w 2 sm y 2nknF 4 ) ■ 


(B37) 


(B33) 
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